NMR investigation of vortex dynamics in Ba(Fe .93Rho.o7)2As 2 superconductor 
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75 As NMR spin-lattice relaxation (1/Ti) and spin-echo decay (I/T2) rate measurements were 
performed in a single crystal of Ba(Feo.93Rho.o7)2As2 superconductor. Below the superconducting 
transition temperature T c , when the magnetic field H is applied along the c axes, a peak in both 
relaxation rates is observed. Remarkably that peak is suppressed for H _L c. Those maxima in 
1/Ti and I/T2 have been ascribed to the flux lines lattice motions and the corresponding correlation 
times and pinning energy barriers have been derived on the basis of an heuristic model. Further 
information on the flux lines motion was derived from the narrowing of 75 As NMR linewidth below 
T c and found to be consistent with that obtained from l/Tb measurements. All the experimental 
results are described in the framework of thermally activated vortices motions. 

PACS numbers: 74.25.nj, 74.25. Wx,74.25.Uv 



I. INTRODUCTION 

The recent discovery of iron-based superconductors! 
was welcomed by the scientific community, because it was 
supposed to answer the still open questions regarding the 
pairing mechanism in high temperature superconductors. 
However the multiplicity of controversial experimental 
results suggests that a unique description of the super- 
conducting properties is far from being reached. Among 
the still debated fundamental topics, e. g. the order pa- 
rameter symmetry^ - — the nanoscopic coexistence of mag- 
netism with superconductivity* 7 -^ the role of antifcrro- 
magnetic spin fluctuations in the pairing mechanism,— ~— 
one fixed point is represented by the study of the flux 
lines lattice (FLL)^ In fact, the study of the magnetic 
field (Tf )-temperature (T) phase diagram of iron-based 
superconductors have immediately attracted lot of inter- 
est owing to their extremely high upper critical fields; 13 ! 14 
which in some cases reach values even larger than those 
of high T c superconductors. 

Most of the theories aiming at describing the FLL 
properties are based on a regular arrangement of 
vortices. 15 However, in real crystals this is far from being 
the case, because crystal defects, such as dislocations or 
inclusions, usually act as pinning centers preventing the 
vortices from having a regular arrangement or from mov- 
ing freely under the action of an electric current. Since 
these dynamics lead to dissipative effects the study of the 
pinning potential is of major importance for the techno- 
logical applications of superconductors. On the other 
hand, the understanding of the different phases devel- 
oping in the magnetic field-temperature phase diagram 
of a superconductor and the modelling of the different 
dynamical regimes give rise to fundamental questions* 16 
A technique which offers the possibility of studying the 
FLL dynamics from a microscopic point of view is nuclear 
magnetic resonance (NMR). In the past years a fruitful 
study was performed in the cuprates 1 ^ - — and showed 
that the linewidth and the spin-lattice relaxation times 
were effective markers of the vortex dynamics. Moreover, 



these two quantities provide complementary information 
since the linewidth narrowing is sensitive to the magnetic 
field fluctuations along the direction of the external field, 
while the spin-lattice relaxation time is sensitive to the 
transverse fluctuations. Furthermore, being the nuclei 
local probes they are sensitive to flux lines excitations 
at all wave- vectors j2£ at variance with macroscopic tech- 
niques, as the AC susceptibility for example, which are 
sensitive just to the long wavelength excitations^ 



Thanks to the works performed in the cuprates we 
know now that in very anisotropic superconductors vor- 
tices can be considered as independent two-dimensional 
isles called "pancakes" which undergo diffusive thermal 
motions! 22 ! 23 Bearing this in mind, and looking at the 
structural similarities between cuprates and pnictides, 
some obvious questions arise: is it still possible to de- 
tect the vortices thermal dynamics in iron-pnictides with 
NMR? What is the vortices structure in the new iron- 
based compounds? Are vortices 2D uncorrelated is- 
lands or rather three-dimensional structures? In order 
to answer at least part of these open questions we per- 
formed an NMR study of the superconducting state of 
Ba(Fei_ x Rh :E )2As2 superconductor with x^ 0.07. We 
measured both the spin-lattice (1/Ti) and spin-echo 
(I/T2) relaxation rates of the 75 As nuclei, together with 
the Knight Shift and the NMR linewidth, at two different 
field intensities (7 T and 3 T) and orientations (H || or 
_L c). The study of these quantities evidences the pres- 
ence of low-frequency dynamics that we interpreted in 
the light of FLL motion and, accordingly, we derived a 
quantitative description of the vortex motions, namely 
the temperature dependence of the correlation time and 
of the pinning potential at different magnetic fields. In 
the present work we will concentrate solely on the super- 
conducting properties, while the discussion of the normal 
state will be presented in a future study. 
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II. TECHNICAL ASPECTS AND 
EXPERIMENTAL RESULTS 

75 As NMR measurements were performed on a flat 
0.8 x 5 x 7 mm 3 parallelepiped shaped crystal of 
Ba(Feo.93Rh .o7)2As2 with the c axis along the shortest 
side. The sample was grown by self-flux method accord- 
ing to a procedure reported in Refll3l. The phase dia- 
gram of Rh-doped compounds shows many similarities 
with that of Co-doped compounds and the maximum es- 
timated transition temperature is about 23 Kii for the 
optimally doped x ~ 0.07 system. For such Rh concen- 
tration both the structural and antiferromagnetic phase 
transitions are suppressed. To provide a first character- 
ization of the crystal we measured the field cooled (FC) 
and zero field cooled (ZFC) magnetization by means of 
a Quantum Design MPMS-XL7 Superconducting Quan- 
tum Interference Device (SQUID) magnetometer. The 
irreversibility line was estimated looking at the temper- 
ature where the ZFC curve departs from the FC one, as 
in ReflH. This is the temperature where the magnetiza- 
tion is sensitive to a change in the dynamics of the FLL. 
On the other hand the detuning of the NMR probed 7 , is a 
higher frequency measurement, which is found to be con- 
sistent with what is observed in static field measurements 
(Fig. P. 



10-1 
8 
6 
4 
2 




o 


T„ (SQUID) 


e 


1/T, peak (NMR) 


* 


TJNMR-probe) 



GO 



18 19 20 21 22 23 24 
T(K) 

FIG. 1: The irreversibility temperature measured with a DC 
SQUID magnetometer (open circle) is compared with that 
derived from the detuning of the NMR probe (blue stars). 
The red circles refer to the temperature of the peaks in 1 /Ti . 
The dotted lines are guide to the eye. 



phase 1/T\T shows a temperature independent behavior, 
as expected for a weakly correlated metal (see the inset 
into Fig. [3]).— By decreasing the temperature below T c 
we observed a well-defined peak in 1/T\ for H || c. The 
peak temperature decreased by increasing the magnetic 
field intensity (see Fig. [3]). Remarkably when H 1 c the 
peak in 1/Ti disappears (see Fig. @|. At lower tempera- 
tures 1/Ti decreases exponentially and it is only weakly 
dependent on the magnetic field orientation. 
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FIG. 2: The recovery curves for three different temperatures 
are shown, for H || c, at 7 T. The blue squares refers to 
the 15 K data, while the pink circles are taken at 17.7 K, in 
correspondence with the peak in 1/Ti, and the black triangles 
refer to 18.7 K. The dotted coloured lines are the best fits 
according to Eq. 1. 
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The NMR measurements were performed by using 
standard radiofrequency pulse sequences. The spin- 
lattice relaxation time T\ was measured by means of a 
saturating recovery pulse sequence at two different mag- 
netic fields H = 7 T and 3 T. The recovery of the nuclear 
magnetization m(t) was found to follow the relatio n 25 ! 26 : 

1 - m(t)/m = 0.1e~ t/Tl + 0.9e~ 6t/Tl (1) 

expected for a nuclear spin / = 3/2 in the case of a mag- 
netic relaxation mechanism (see Fig. [5]). In the normal 



FIG. 3: The spin-lattice relaxation rate, measured at 7 T 
(open squares) and 3 T (blue circles), for H || c is reported. 
The inset shows the 1/TiT data at 7 T both in the supercon- 
ducting and normal phase. The arrows show the temperature 
of the detuning of the NMR probe at the two fields: the blue 
arrow stands for 3 T and the black arrow for 7 T. 

The transverse relaxation time was measured by 
recording the decay of the echo after a7r/2 — r — ir pulse 
sequence as a function of the delay r. Since the functional 
form of the decay changes with temperature (Fig. [3]) , as 
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FIG. 4: The spin-lattice relaxation rate, measured at 7 T, 
in H || c geometry (black diamonds) and H _L c geometry 
(blue circles) is shown. A neat difference for the two field 
orientations is found in the 16-19 K range. Data, in H 1 c 
geometry, have been normalized by a value 1.55 to match the 
value of 1/Ti for H || c, at T c thus revealing an anisotropy of 
the hyperfine tensor. 



will be discussed subsequently, in order to compare the 
data over the full temperature range we defined Ti as the 
time where the echo amplitude decreases by 1/e. 



spin echo decay rate is found to reach the value derived 
from nuclear dipolar lattice sums. 

Similarly to what was observed for 1/Ti, also I/T2 peak 
gets significantly reduced for H _L c (Fig. [5]). 




FIG. 6: The figure shows the spin echo decay rate measured at 
7 T. A peak around 12 K is found for H || c (black diamonds) 
while it strongly decreases, for H _L c (blue squares), and 
shifts towards higher temperatures. The red arrow indicates 
the ab-initio value for I/T2 given by the dipolar sums. The 
inset shows the spin echo decay rate for H || c up to room 
temperature. 




6 2 4 6 8 10 

2t (ms) 



FIG. 5: The figure shows two echo decays as a function of 2r, 
below Tdet '■ the black squares refer to 24 K and the blue circles 
refer to 15 K, below the superconducting transition where the 
FLL is still dynamic. From the figure one can notice that the 
echo decay functional form changes while the temperature 
decreases. The red curves are the best fits according to Eqs. 
10 and 11. 

In the normal phase 1 /T2 shows an activated tempera- 
ture dependence whose origin will be discussed elsewhere. 
Below T c we observed a marked increase in I/T2 giving 
rise to a peak around 12-13 K, for H || c. We note that 
Oh et al. found a similar behaviour in their data referred 
to the 7.4% Co-doped single crystal.— Nevertheless we 
note that the compound is different though 1/T 2 data 
are quite similar: they observed a peak in I/T2 around 
15 K while we observed it around 12-13 K. For T — > the 



The NMR spectrum was determined from the Fourier 
transform of half of the 75 As echo signal, while below 
T ~ 13 K, when the line became too broad, the spectrum 
was derived by sweeping the irradiation frequency. The 
full width at half maximum (FWHM) was determined 
by a Gaussian fit. In the normal state the linewidth in- 
creased on cooling following a Curie- Weiss trend (Fig. [7]) , 
probably due to the presence of paramagnetic impurities. 
The impurities cause the appearance of a staggered mag- 
netization and a broadening of the NMR line. On the 
other hand, the average magnetic field is only weakly af- 
fected, so we do not expect an extra-contribution to the 
shifts After subtracting this impurity-dependent con- 
tribution Ai>np from the raw data by using the relation 

Ai/(T) ~ y/Av(T)* aw -A^ NP (2) 

we observed that below T c an extra-broadening induced 
by the presence of the flux lines lattice appears (Fig. [7]). 
The impurity-dependent contribution was found to be 
well described by the Curie- Weiss relation, for both the 
sample orientations: 

A VNP (T) = ^— e +A (3) 

By assuming the value of 9 = —60 K, as found by the 
fitting procedure we obtained the following results: for 
the H || c case the fit gave C = 1319 ± 118 kHz K and 
A = 20.8 ± 1 kHz, while for the perpendicular geometry 
the fit gave C = 1264 ± 50 kHz K and A = 23.4 ± 0.1 
kHz. 
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FIG. 7: The figure at the top shows the FWHM at 7 T, for 
H || c, after the Curie- Weiss correction (see Eq. [2]). An effect 
of narrowing occurs just beiow T c : the experimental data 
deviate from the London two- fluid model (blue line). The 
figure at the bottom shows the FWHM at 7 T, for H _L c, 
after the Curie- Weiss correction. Here the extra-broadening 
can be fitted by the London two-fluid model, up to T c . In the 
insets the Curie- Weiss behavior observed for the raw data is 
shown. 



It has to be noticed that the superconducting state 
affects not only the 75 As NMR linewidth but also the 
NMR shift. Above T c , in the normal phase, the NMR 
shift shows an activated behavior, as observed also for 
the Co-doped BaFe2As2.— The experimental data (Fig. 
|5]) can be fit with an activated Arrhcnius law: y = A + 
5exp(-T>/T), yielding A = 0.26 %, B = 0.071 % and 
D = 225 ±22 K, for H_JI c, in good agreement with 
the values found in Refl28l. Below T c the shift starts 
to decrease as expected for a singlet state pairing^ In 
the superconducting phase the NMR shift K(T) can be 
assumed to result from three contributions: 



K{T) = K spm (T) + K FL {T) + K TI 



(4) 



shift, orbital terms, etc...)4i Owing to the line broad- 
ening and to the reduction in the radio-frequency pene- 
tration depth, the accuracy in the estimate of K(T) de- 
creases below T c and does not allow us to draw convincing 
conclusions on the symmetry of the order parameter. 
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where K sp i n (T) is the spin-dependent part, which van- 
ishes for T — > 0, Kfl{T) is the diamagnetic correction 
due to the vortex lattice, and the last term contains 
all the temperature independent contributions (chemical 
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FIG. 8: The figure shows the Knight Shift at 7 T, for H || c 
with an Arrhenius-like fitting curve (solid line) in the normal 
phase. 

Taking into account the quadrupolar shift of the cen- 
tral line, in the transverse geometry, we estimated a 
quadrupolar frequency ^q(100 K) ~ 1.5 MHz, very close 
to the one found in the parent compound.— 



III. DISCUSSION 

As previously mentioned, here we will not discuss the 
normal state properties but rather we shall concentrate 
on the superconducting phase. Let us first consider 
the behaviour of the spin-lattice relaxation rate which 
is characterized by a well defined peak for H || c be- 
low T c . In passing, we note that in Co-optimally doped 
compound^ no peak was observed in 1/Ti, below T c . 
On the other hand Laplace et air- in a 6% Co-doped 
BaFe2As2, found a peak in 1/T\ just below T c and an 
enhancement in \jT\T at higher temperature due to spin 
density wave correlations. The peak we found below 
T c is not expected to be a Hebel-Slichter peak^ since 
the majority of the experimental and theoretical results 
point towards an extended s -wave pairing)2iir— ^ where 
that feature is expected to be absent. Furthermore if 
it was a coherence peak the data would be described, 
below T c , by 1/Ti ~ e~ A / T with A the superconduct- 
ing gap. By fitting the data one obtains A ~ 200 K 
>> 3.5fcsT c , the value expected for the superconduct- 
ing order parameter^ Finally, the striking suppression 
of the peak for H 1 c can hardly be reconciled with the 
small anisotropy of the electron spin susceptibility found 
in those materials. Hence that maximum in 1/Ti just 
below T c should not be associated with the electron spin 
dynamics but, given the similarities with the behavior 
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found in HgBa 2 Cu04 +1 j2 6 . and YBa 2 Cu40g2 7 . cuprates, 
it is tempting to associate 1 /Ti peak to the FLL dynam- 
ics. 

In order to analyze the experimental results one 
can start from the basic modelling of FLL in strongly 
anisotropic superconductors^ the vortices enter the 
sample in form of quasi-two dimensional pancakes, ly- 
ing in the FeAs planes. Owing to the thermal excitations 
they move out of their equilibrium positions by means 
of random motions, which can be hindered by the pin- 
ning centers. Differently from cuprates, which exhibit 
a very high anisotropic ratio 7 = £ a b/£ c , with ^ a b and 
£ c the in-plane and out of plane coherence lengths, the 
Bal22 superconductors show 7 ~2-4 varying with tem- 
perature^ This suggests to describe the flux lines not as 
completely uncorrelated pancakes, but rather as a stack 
of correlated islands. Still, since the estimated correla- 
tion length £ c is of the order of the inter-layer distance 
Sj22, namely 2£ c ~ s ~ 6 A , FeAs planes can be considered 
as weakly coupled superconducting layers. Accordingly, 
when H 1 c the flux lines are preferentially trapped be- 
tween the planes and the FeAs plane boundaries act as 
pinning centers, a well known effect in layered supercon- 
ductors. These intrinsic pinning centers hinder the dy- 
namics and yield the observed suppression in the 1/Ti 
peak for H 1 c. 

In order to understand the shift of the 1 /Ti peak upon 
increasing H we first recall that 1/Ti probes the spec- 
tral density J(u>l) at the nuclear Larmor frequency u)l,, 
namely 



1 7 2 f 

Y 1 =TJ < M*)M0) > e 



iujLt 



(It 



(5) 



with hp the magnetic field component perpendicular 
to H and 7 = 2tt x 7.292 x 10 6 rad/T the gyromagnetic 
ratio of the 75 As nucleus. Then the field dependence of 
the peak in 1/Ti can be qualitatively understood by con- 
sidering that at the peak temperature the characteristic 
frequency for FLL motions is close to ujl- When the mag- 
netic field increases, T c decreases and so does Ti rr , hence 
the FLL dynamics remain fast over a broader tempera- 
ture range and the maximum in 1 /Ti is observed at lower 
temperature (Fig. 1). It is noticed that the peak in the 
spin-lattice relaxation rate appears just below the irre- 
versibility temperature, in contrast with what was found 
in the cuprates, where it is well below the irreversibil- 
ity line, suggesting a higher FLL mobility in these latter 
compounds j^ 7 - 

To give a quantitative description of the peak we 
started from the equation [5j Let us first assume that the 
vortex fluctuations are basically two-dimensional (2D), 
that take place in a spatial range smaller than the 



inter- vortex distance 4 ^ \ e = y2/y/3^^o/H (for a tri- 
angular FLL) and that they move by Brownian mo- 
tion o 36 ' 37 described by a diffusive-like correlation func- 
tion g\(t) = exp(—Dj_q^t) 7 being the diffusion con- 
stant of the motion taking place in the ab plane. Then 



T c {q±) = 1/D±_q\ plays the role of a q-dependent cor- 
relation time for the collective vortex motions. By sum- 
ming over all collective in-plane excitations up to a cut 
off wave- vector q m = (l// e )(8Tr 3 /3) 1 / 4 Suh et al. 36 found 
the spectral density 



J (ujl) = r m In 



(6) 



where the average correlation time is r m — 1/ D±qf n . For 
the temperature dependence of r m it is reasonable to as- 
sume an activated form r m (T) = r exp(U/T), where U is 
an average pinning energy barrier and To stands for the 
correlation time in the infinite temperature limit. Ac- 
cordingly an activated temperature dependence of the 
spectral density at the Larmor frequency and then of 
1/Ti are observed for T — > 0. The best fits of the data 
according to this 2D vortex model are reported in Fig. 
8. It is noticed that the fit is not fully satisfactory. 

On the other hand, as previously pointed out, the low 
anisotropy of BaFe 2 As2 compounds suggests that signif- 
icant vortex correlations along the c axes are present in 
Ba(Feo.93Rho.o7)2As2- Thus, the flux lines have to be 
considered as stationary waves oscillating in between the 
pinning centers. In order to take into account this effect 
we introduced empirically a modulation in the amplitude 
of the correlation function characterized by a wavelength 
A which has an upper bound given by A c , the London 
penetration depth along the c axes. Then one can write 
<?2(£) = exp(— D±q^t) cos (z/\). Now if we recall the 
form of the longitudinal field correlation function 37 



< h p (0)h p (t) > 
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and taking the root mean square amplitude of the vortex 
core fluctuation with respect to equilibrium positio n 41 ! 42 
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Eq. [7] can be written: 
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where the temperature dependence is evident. Taking the 
values for the London penetration depth reported in the 
literature^ 3 - the coherence lengths derived from the H c2 
measurements, and their temperature dependence ac- 
cording to the two-fluid model, we were able to reproduce 
fairly well the temperature dependence of 1/Ti, which 
indicates that indeed a " 3D-correlated- vortices" model 
is more appropriate to describe Ba(Fei_ 2: Rh 2: )2As2 su- 
perconductors. The best fits of the experimental results 
(Fig. 1]) gives a value of U = 322 ± 66 K for H = 7 T 
and U = 470 ± 5 K for H = 3 T. These values are similar 
in magnitude to those found in YBCO-124^ 7 - nonethe- 
less the quality of the fitting procedure suggests that the 
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vortices develop a three-dimensional correlation. Before 
concluding this part we estimate the root-mean-square 
amplitude of the transverse field fluctuations < h 2 e > 
which represents the ripple of the magnetic field profile 
modulated by the flux lines dynamics. Infact given the 
Eqs. [5] and 13 1/Ti can be written in this new form 
1/Ti = (j 2 /2) < h 2 > J(lol) from which we obtained 
h e ~ 30 - 40 Gauss at 7 T and - 20 Gauss, at 3 T. We 
point out that these values are close to the low temper- 
ature NMR full width at half maximum, as it has to be 
expected. 
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FIG. 9: The figure at the top shows the spin-lattice relax- 
ation rate at 7 T for B || c, while the one at the bottom 
shows the spin-lattice relaxation rate at 3 T in the same ge- 
ometry. The fitting curves are given by the 2d - uncorrelated 
pancakes model, deriving from the correlation function gi(t) 
(dash-dotted line) and the correlated vortices model, deriv- 
ing from the correlation function 512 (t) (solid line). In both 
cases the second model shows the best agreement with the 
experimental data. 

While the spin-lattice relaxation rate has been con- 
sidered one of the most valuable microscopic probes of 
the FLL motion, not so much effort has been devoted to 
the analysis of the spin echo decay time, mainly because 
its interpretation is not always straightforward , 44 i 45 As 
it has been already pointed out in the superconduct- 
ing state I/T2 shows a neat anisotropy: the peak found 



for H || c is significantly reduced and shifted in the 
H _L c configuration. Moreover, we point out that at 
low-temperature I/T2 reaches the value expected from 
the Van Vleck lattice sums i 28 ' 46 In fact, at low tempera- 
tures all dynamics are frozen and the only process giving 
rise to the echo decay is the nuclear dipole-dipole inter- 
action. The peak observed around 10 K cannot be due 
to a time-dependent modulation of that interaction since 
we do not expect such an anisotropic behavior, in that 
case. On the other hand, given the similarity with I/Ti 
peak anisotropy, it is likely that also I/T2 peak arises 
from a low frequency vortex dynamics. Nevertheless it 
should be emphasized that while I /T\ measurements are 
sensitive to the fluctuations of the transverse components 
of the magnetic field, I/T2 is sensitive to the longitudi- 
nal ones. In particular, it should be noticed that when 
the vortices are strongly correlated along the c axes, the 
flux lines move rigidly and do not affect significantly the 
transverse field components, while they do change the 
longitudinal ones 42 Hence the information derived from 
those two types of measurements can be complementary. 

In order to analyze the temperature dependence of 
I/T2 we need an analytical expression for the spin echo 
decay. In principle we should start from a relation simi- 
lar to Eq. [7J nevertheless here, for the sake of simplicity, 
we assumed an exponential correlation function for the 
longitudinal field fluctuations 

< hi{Q)hi{t) >=< hf > e~ t/TL (9) 

characterized by an average correlation time T£. Corre- 
spondingly we have written the decay of the echo ampli- 
tude as a function of the delay r between the ir/2 and tt 
pulses in the echo sequence as^£ 

M(2t) = M e~ 2r/T ^" x M 2 (2r) (10) 

M 2 (2t) = e -7 2 <'j?>T£[2r/r I .+4o X p(-2r/r I .)-o X p(-2T/T L j) r j3j 

where the first Gaussian term accounts for the nuclear 
dipole-dipole contribution, while the second term de- 
scribes the low-frequency vortex motions. By fitting the 
data below T c we were able to derive the temperature de- 
pendence of the longitudinal correlation time (Fig. ITU1) . 
By decreasing the temperature the FLL motion is sup- 
posed to go through different motional regimes: from 
the fast motions (tl *C T2) up to the very slow motions 
(t£ ^ T2), where the correlation time is so long that we 
can consider the FLL to be frozen in the solid state. If 
we fit the data above II K, where the peak in 1/Ti is ob- 
served, we notice that tl follows an activated behavior 
characterized by an activation barrier Ul ~ 50 K much 
lower than the one derived from 1/Ti (see. Fig fTOj) . Here 
we refer just to the fast motion limit, since the fitting pro- 
cedure is not so much accurate at the low temperatures, 
because of the reduced signal to noise ratio. 

The correlation time derived from the spin-echo decay 
rate can be suitably compared to that derived from the 
motional narrowing of the NMR line. The latter can 
be derived, for the H || c case, following the standard 
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approach reported in Ref. |4 
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In the fast motions regime, 

1/2 

tl « 1, with (A^) 2 the square 



namely 27t(Az;r) 2 
root of the rigid lattice second moment, one has 



Av ~ T L 



(A./ fl ) 
2vr 



(12) 



By fitting the correlation time T£ with an Arrhenius's 
law we extracted a pinning energy barrier Ul — 48 ± 3 
K, which is consistent with the one derived from the spin 
echo decay measurements. This is not surprising since 
both Ti and Av probe the longitudinal component of 
the local field fluctuation. 

In order to understand why the energy barrier probed 
by I/T2 and from the motional narrowing are smaller 
than the one derived from 1/Ti measurements it should 
be pointed out that the oscillations at wavevector qu — > 
do contribute to the longitudinal field fluctuations but 
only weakly to the transverse field excitations which are 
relevant in 1/Ti, at variance with the 311 — > 1/s modes, 
which contribute significantly to 1/Ti. Hence, our find- 
ings indicate that the energy cost to activate a certain 
collective mode increases with increasing qu, where the || 
subscript refers to the wavevector component parallel to 
the magnetic field. 
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FIG. 10: The temperature dependence of the correlation time 
tl gives by echo decay time (black squares) is compared with 
the one derived by the linewidth analysis (blue circles) in the 
assumption of fast motions (see Eq. [12}. The red curves are 
the fitting of the correlation time, according to an activated 
law. 



time gets longer than a few ms. The Gaussian lineshape, 
instead of the asymmetric one expected for a perfect tri- 
angular lattice, indicates the presence of lattice distor- 
tions induced by randomly distributed pinning centers. 
In this scenario one can only make an estimate of the 
London penetration depth which can be compared with 
the one derived by transverse fiSK^ transport 43 and the 
tunnel diode resonator measurements^ on a similar 7.4% 
Co-doped BaFe2As2 single crystal. These authors report 
X a b values between 200 and 217 nm. Following Ref. [5l| 
we extracted: 



'2. 36^07^ 



(13) 



where \fk = 0.04324 depends on the lattice geometry and 
on the magnitude of the applied field. Taking Av(T — >• 
0) ~ 30 kHz, for H || c we extracted A a6 (0) - 226 ± 9 
nm, in agreement with the former results. 



IV. CONCLUSIONS 

This work aimed at studying the thermally activated 
vortices motion, by means of 75 As microscopic probe. We 
found that in the NMR relaxation rates and in the NMR 
linewidth there is evidence of such motions, which is well 
supported by the remarkably anisotropic behaviour of 
those quantities. We were able to follow the dynamics of 
the vortices by measuring NMR quantities which are sen- 
sitive to motions at different time scales: the peak in the 
spin-lattice relaxation time is found when the correlation 
time is about 10" 7 - 10~ 8 s, while the 1/T 2 maximum oc- 
curs at a slightly lower temperature, when the correlation 
time is comparable to I/T2, i.e. few ms. In the tempera- 
ture window between those peaks the motions are effec- 
tive and yield the motional narrowing of the NMR line. 
To observe a line narrowing the correlation times must 
be smaller than the inverse of the rigid lattice linewidth 
~ 10~ 4 s. Furthermore we pointed out that the tem- 
perature dependence of 1/Ti around the peak suggests 
that the flux lines are formed by strongly coupled vor- 
tices rather than nearly independent pancakes diffusing 
in 2D, as it was found in the cuprates. 
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